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The scaling of island and monomer density, capture zone distributions (CZDs), and island size 
distributions (ISDs) in reversible submonolayer growth was studied using the Clarke- Vvedensky 
model. An approach based on rate-equation results for irreversible aggregation (IA) models is 
extended to predict several scaling regimes in square and triangular lattices, in agreement with 
simulation results. Consistently with previous works, a regime / with fractal islands is observed at 
low temperatures, corresponding to IA with critical island size i = 1, and a crossover to a second 
regime appears as the temperature is increased to eR 2 ^ 3 ~ 1, where e is the single bond detachment 
probability and R is the diffusion-to-deposition ratio. In the square (triangular) lattice, a regime 
with scaling similar to IA with i = 3 (i = 2) is observed after that crossover. In the triangular lattice, 
a subsequent crossover to an IA regime with i — 3 is observed, which is explained by the recurrence 
properties of random walks in two dimensional lattices, which is beyond the mean-field approaches. 
At high temperatures, a crossover to a fully reversible regime is observed, characterized by a large 
density of small islands, a small density of very large islands, and total island and monomer densities 
increasing with temperature, in contrast to IA models. CZDs and ISDs with Gaussian right tails 
appear in all regimes for R ~ 10 7 or larger, including the fully reversible regime, where the CZDs 
are bimodal. This shows that the Pimpinelli-Einstein (PE) theory for IA also explains the main 
mechanisms for the large islands to compete for free adatom aggregation in the reversible model, 
and may be the reason for the successful application of that theory to a variety of materials and 
growth conditions. 

PACS numbers: 68.43. Hn, 05.40.-a, 68.35.Fx , 81.15.Aa 



I. INTRODUCTION 



The morphology and physical properties of thin films 
and multilayers are strongly related to the initial stages 
of their formationir— . It motivated a large number of 
works on modeling the submonolayer growth regime^—, 
in which a single incomplete layer of adatoms is being 
formed. These models usually consider a set of funda- 
mental processes (diffusion, aggregation, etc) in order 
to find activation energies and related quantities for a 
particular deposition process or with the aim of investi- 
gating universal features valid for several materials and 
techniques. 

The simplest models consider irreversible aggregation 
(IA) of atoms to islands of size larger than a critical value 
i. Rate-equation (RE) approaches predict the scaling 
of island and monomer densities (iVj a {, Ni)££. for any 
i and show good agreement with simulation data, usu- 
ally performed for i ranging between 1 and 3*. The 
shapes of the island size distributions (ISDs) were pro- 
posed by different approximation schemes^— and are still 
subject of debate. A recent advance by Pimpinelli and 
Einstein (PE)^S proposed that the capture zone distribu- 
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tions (CZDs) are described by the Wigner surmise (WS) 
from random matrix theoryii [the capture zone (CZ) 
of an island is defined as the area in which a diffusing 
adatom is more likely to attach to that island than to 
any other one]. After an initial controversy^— , that 
proposal was supported by simulation results for vari- 
ous island shapes after suitable rescalingJ^, under condi- 
tions of high temperatures and low coverages. IA models 
are expected to model real systems in temperatures suffi- 
ciently low to neglect atom detachment from islands, but 
this condition is not obeyed in many real systems, even 
at room temperature. For this reason, suitable choices of 
the size i are frequently combined with internal restruc- 
turing mechanisms to reconcile experimental data and IA 
prediction o 7 ' 16 " — . 

On the other hand, a smaller number of works 
have addressed scaling properties of models with re- 
versible aggregation (RA)22r— . Works on square lattices, 
which simulate metal (100) homoepitaxy, have shown a 
crossover between IA scaling regimes with i — 1 and i = 
3 21 i 23 . Partially reversible models allowing only single- 
bond detachment were studied using simulations and im- 
proved RE approaches and confirmed the crossover scal- 
ing between those regime o 27 ' 28 . The work on triangu- 
lar lattices, which is related to (111) epitaxy, shows a 
crossover between regimes with i = 1 and i = 2 27 i 28 . 

In this paper, we study the scaling properties of island 
density, CZDs, and ISDs in a model of fully reversible 
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island growth, using numerical simulations in square and 
triangular lattices and a scaling approach based on RE 
results for IA models. We consider a canonical bond- 
counting model of Clarke- Vvedensky22 type, in which the 
activation energy for an adatom hop to a neighboring site 
depends only on the number of nearest neighbors at its 
initial position, with the temperature T kept fixed dur- 
ing the growth. It respects detailed balance conditions, 
in contrast to IA. For a fixed coverage 9 and increasing 
deposition temperature, the islands evolve from a frac- 
tal shape (typical of IA with no island restructuring) to 
a compact shape dominated by processes of adatom de- 
tachment and reattachment, and the island density con- 
tinuously decrease. The previously studied crossovers of 
IA regimes are reproduced, but an additional crossover 
from i — 2 to i = 3 in triangular lattices is found. At 
high temperatures, a nontrivial regime where the island 
density increases with temperature is observed, with a 
large density of small clusters coexisting with large com- 
pact islands. The Gaussian right tail of CZDs predicted 
by PE theory is observed in all regimes for small cover- 
ages and sufficiently high temperatures, even when the 
CZDs are bimodal. 

The rest of this work is organized as follows. In Sec. [Ill 
we present the RA model and summarize previous theo- 
retical results for IA and RA. In Sec. IIII1 we present the 
simulation results and a scaling approach in square lat- 
tices. In Sec. IIV1 that discussion is extended to triangu- 
lar lattices. Sec. IVl summarizes our results and presents 
our conclusions. 



II. BASIC DEFINITIONS AND THEORETICAL 
APPROACHES 

In the RA model, the adsorbed atoms can occupy sites 
of a square or triangular lattice, with the maximum of 
one adatom per site. There is a random flux of F atoms 
per site per unit time. Adsorption is allowed only if the 
site of incidence is empty, so that a single layer is formed 
on the substrate lattice. This condition is reasonable for 
small coverages. The hopping rate of an adatom at a 
given site f is given by D = Dq cxp (— E act /kBT) (pro- 
portional to its surface diffusion coefficient), with the ac- 
tivation energy E act = E s + uEn , where n is the number 
of adatoms in nearest neighbor (NN) sites of r, E s is an 
energy barrier for surface diffusion of free adatoms, and 
En is the bond energy for each NN. 

Simulations were performed in very large lattices, typ- 
ically of lateral size L = 1024, with at least 100 realiza- 
tions up to coverage 6 — 0.4 for each parameter set. 

The parameters F = 0.1s -1 , D = 10 13 s _1 , and E s = 
1.3eV were kept fixed in all simulations. The last value 
is characteristic of Fe/Fe(100) epitaxy 20 . Broad ranges 
of temperature are studied, thus the conclusions can be 
extended to other systems by considering the diffusion- 
to-deposition ratio of free adatoms R = D/F as the main 
model parameter. Due to the large frequency Dq, one 



typically has R ^> 1, even at low temperatures. 

Several values of the bond energy En are analyzed, and 
scaling approaches are facilitated by using the detach- 
ment probability (binding factor) e = cxp (— En /hsT). 
Since En > 0.3ey is considered throughout this work 
and the largest experimentally relevant temperatures do 
not exceed fcgT ~ O.leV, the condition £ < 1 always 
applies. It is also important to observe that both R and 
e increase with temperature. 

In the I A models, islands with a number of atoms 
larger than the critical size i are stable. Moreover, islands 
of size i are assumed to have a binding energy Ei < 0. 
During the deposition process, after a transient regime, 
islands grow by capture of diffusing adatoms, with negli- 
gible formation of new islands. In this steady state, rate- 
equation theory^ predicts that the stable island density 
scales as 

N i8l ~ e 1 '^ exp {-pEi/ (i + 2)]R-* (1) 
and the monomer (free adatom) density scales as 

Nt ~ 0- 1 /( J + 2 ) exp \pEi/ (i + 2)]RX- 1 , (2) 

where 

X = if (i + 2) . (3) 

These results assume that the mean capture number of 
stable islands is independent of the coverage. 

Information on the growth dynamics can also be ex- 
tracted from the probability densities of island size s, 
Q(s), and of CZ area x, P(x) (ISD and CZD, respec- 
tively). The ISD follows the scaling form 



Q00 = yVfrv 



(4) 



where / is a scaling function, and an equivalent scal- 
ing form applies to P(x). Alternatively, scaling with the 

2 1/2 

variance a x = (x — (x)) may be used as 



P(x) 



1 fx— (x) 

—g 



(5) 



This procedure was previously used with IA models^ and 
is helpful in other problems, such as scaling of roughness 
distributions of thin film growth models^. 

The best know theory to predict the shape of CZD in 
IA is that of PE, which proposes the CZD to be described 
by the WS 



Pp{z) = apz^exp (-bpz 2 ), 



(6) 



where z = xj (x), /3 = | (i+ 1), d is the substrate di- 
mension (d = 2 in the present work) and the parameters 
ag and bp are determined by normalization conditions. 
This theory follows from the phenomenological argument 
that the CZD can be extracted from a Langevin equa- 
tion representing the competition of neighboring islands 
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FIG. 1: Square lattice islands for (a) T = 650K, (b) T = 900K, and (c) T = 1200if. The panels are 100 x 100, 200 x 200, and 
300 x 300 cuts of systems of size L = 1024, respectively. The coverages are 8 = 0.05, 0.10, 0.20, and 0.30 from left to right. 



for adatom aggregation. Several experimental works 
have already shown agreement of CZDs with the WS, 
such as growth of para-sexiphenyl islands^ 3 -?, Cu depo- 
sition with impurities 3 - 3 -, pentacene island growth with 
impurities^, InAs quantum dot growth on GaAa^, and 
Cqo deposition on SiC>2 films2&. 

In Ref. 15 , a scaling approach was used to predict the 
decay of the right tail of the ISD from the Gaussian tail 
of the CZD (Eq. |6]) and the island shape (point, fractal 
or square). 



III. SQUARE LATTICE 

Figs, la-c illustrate the island evolution for Em — 
OAeV and temperatures representative of different scal- 
ing regimes. Fig. 2a shows the temperature evolution 
(parametrized by R) of the island density and monomer 
density for coverage 8 = 0.1 and several values of En- 

For very low temperatures (R < 10 3 ), both densities 
are nearly temperature-independent, and the system be- 
haves as in random sequential adsorption (RSA) without 



diffusion 3 - , as discussed in Refits. 

The subsequent low temperature regime, hereafter 
called regime 7, is characterized by the decrease of Nui 
and Ni with R, as shown in Fig. 2a. In this regime, 
the probability of detachment of an aggregated adatom 
from an island is negligible during the time interval nec- 
essary for the aggregation of a new adatom to one of its 
neighbors. Thus, the average detachment time of singly 
bonded atoms, 

n~l/0De), (7) 

is much larger than the time for aggregation of a new 
atom at a typical island. In a mean-field approach (RE), 
the latter is given by 

r ag ~ 1/ (DNt) . (8) 

Regime I is equivalent to IA with i = 1, and using Eq. 
([1]) for fixed coverage one has 

N lsl ~ iT 1 / 3 , Ni - R-' 2/3 (I) (9) 
This scaling is highlighted in Fig. 2a. 
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FIG. 2: (Color online) Scaled island densities for 9 — 0.10 
and several binding energies in the square lattice. Dashed 
lines indicate local slopes. The insets show scaled monomer 
densities. 



Despite the effective IA scaling of regime i, the CZDs 
usually do not match the WS with i = 1 (P 2 ) for the 
values of En and temperatures considered here. Fig. 
2a shows that regime I typically occurs with i? < 10 7 . 
This range is intermediate between low adatom mobility 
(RSA behavior) and high adatom mobility (PE theory), 
thus CZDs show a crossover scaling, as shown in Ref. 38 . 
Similar crossover is observed in the ISDs. For instance, 
considering Ejy = 1.0, for i? < 10 5 we observe simple 



exponential tails of CZDs (RSA scaling) and for R ~ 10 7 
they become Gaussian (PE theory), but still not collaps- 
ing with the WS. 

A crossover to a second regime is observed when 
t\ matches r ag [Eqs. (JTJ) and ©], as explained in 
Refsj 21 i 23 ' 27 (see also Sec. 8.5 of Ref4). A crossover 
variable is defined as 



Y x = eR 2/3 , 



(10) 



with Y\ <C 1 in regime i. Regime II begins with Y\ ~ 10. 
From Eqs. ^ and (fTU|) . the crossover scaling for island 
density and monomer density is 



N„ 



is I 



(11) 



where F\ and G\ are scaling functions. Regime II is 
effectively characterized by i = 3, since detachment of 
doubly bonded atoms of square islands (size s = 4) occurs 
in much longer timescales. Thus, for Y\ > 10 but not very 
large, Eqs. (fT]). (|2|), and (|TT]) give 

N lsl ~ e-V 5 R-V\ Ni ~ e 2 / 5 i?- 2 / 5 (// initiol ). 

(12) 

This result is consistent with a binding energy E, = 2En 
for the critical island with i = 3 adatoms. Fig. 2b illus- 
trates the crossover scaling with the collapse of data in 
regimes I and 77 for various E jy. 

For En = 1.0, Fig. 2a shows that the slope of the 
Ni S i x R plot slowly increases in regime II. This occurs 
because the binding energy Ei becomes negligible at the 
end of this regime, so that Ni S i and Ni scale tend to scale 
only with R, as 



iT 3 / 5 , A^-iT 2 / 5 {Ilfinal) 



(13) 



[the same power laws of Eq. (fT2|) excluding the e- 
dependence]. The same evolution of slopes is observed 
in all effective IA regimes observed in the RA model (see 
Sec. iTVTt . 

As the temperature increases, a crossover to a new scal- 
ing regime is expected by two reasons: rapid detachment 
of doubly bonded atoms from islands with s > 4, which 
occurs in a characteristic time 



T-2 



1/ (De 2 ) , 



(14) 



or slow atom deposition in the CZ, which occurs in a 
characteristic time 



Tdep ~ Nisl/F. 



(15) 



In the steady state of IA, Tdep [Eq. (|T5|) ] is of the same 
order of r ag [Eq. ©]. Thus, regime // is expected to end 
as t 2 [Eq. flip] matches T ag , i. e., when ei? 1 / 5 ~ 1. This 
defines a second crossover variable as 



Y 2 = ei? 1 / 5 . 
Crossover scaling follows from Eq. (fT3")l as 

N M ~ e 3 F 2 (Y 2 ) , N± ~ e 2 G 2 (Y 2 



(16) 



(17) 
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FIG. 3: (Color online) Scaled (a) CZDs and (b) ISDs for the 
energy En = 0.65 and temperatures T = 950A (7? ~ 10 7 ) 
and T = 1150if (R ~ 2 x 10 8 ). The insets show the same 
scaled distributions with squared abscissas. The dashed lines 
highlight the Gaussian decay of the right tails. 

Fig. 2c illustrates the crossover scaling for various En, 
with an excellent collapse of N is i data. Some deviations 
appear for 7Vx, probably due to scaling corrections. 

After this crossover, a fully reversible regime III is 
attained, completely different from IA. Fig. 2a shows 
that the island density and the monomer density increase 
with temperature in this regime, in striking contrast with 
regimes 7 and 77 and all IA models. Fig. lc shows that 
it is characterized by a very large density of small is- 
lands, mainly isolated adatoms, and a small density of 
very large islands. 

In regime 777, the density of small islands is much larger 
than the density of large islands. This is an expected evo- 
lution from regime 77 because N\ decreases with temper- 
ature slower than Ni S i [Eq. (|12l) ]. This trend is expected 
for any regime with i > 2, which gives \ > 1/2, and will 
also be shown in the triangular lattice. 

The high-temperature scaling in regime 777 is only e- 
dependent, with 

A^~e 3 , Aq~e 2 (777). (18) 

Indeed, since movement of free adatoms is very fast, only 



FIG. 4: (Color online) (a) Scaled ISDs for E N = 0.30, T = 
950 A' (7? « 10 7 ) and several coverages. The inset shows the 
same data in a log-log scale and the line has slope —4. (b) 
Scaled CZDs for the same parameters of (a). The inset shows 
the same distributions with squared abscissas. The dashed 
lines highlight the Gaussian decays of the right tails. 



the attachment-detachment dynamics is important to de- 
termine the densities. 

This steady state regime can be predicted by noting 
that double bonding is the minimal possible bonding for 
an island in the square lattice, i. e. at least four atoms 
of a large island have two bonds. Equating aggregation 
and detachment rates obtained from Eqs. JHJ and (|14[) . 
we get the above relation for N\. Since regime 777 takes 
place as the evolution of an effective i = 3 IA regime 
(also related to lattice structure), where Ni s i ~ N^ 2 , we 
obtain the above relation for Ni S i . 

The CZDs in regime 77 are not well fitted by the WS 
with i = 3 (7-4 in Eq. [5]), probably because this is not 
a true IA process, as shown in Fig. 3a. However, as 
the temperature increases, the CZD gets closer to that 
WS. This occurs considering scaling with the average [Eq. 
dU] or with the variance [Eq. ((5)]. On the other hand, 
Fig. 3a shows that the tails of the CZDs are Gaussian, 
similarly to IA models^. This trend fails only for very 
small 7?jv because regime 77 is attained with small 7? and 
crossover from RSA is present. 
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The ISD has features similar to those of IA models 
of compact islands (e. g. square islands), as shown in 
Fig. 3b: the left tail is high, indicating the presence of 
a large number of small islands, and the right tail shows 
a Gaussian decay. As the coverage increases, a crossover 
to simple exponential decay is observed in the CZDs and 
ISDs, also consistent with the IA models 38 . 

Fig. 4a shows ISDs in regime III, with left tails much 
higher than the peaks, confirming the presence of a large 
number of free adatoms and small islands. The left tail 
shows even-odd oscillations characteristic of loose-packed 
lattices (inset of Fig. 4a). If the small islands were 
in near-equilibrium conditions, the Walton relation^i 3 ^ 
Q(s) » Ni s would apply and an exponential decay of 
the left tail would be observed. However, this is not 
the case: the inset of Fig. 4a shows an approximately 
power-law decay, with large exponents that depend on 
the coverage, suggesting that more complex mechanisms 
govern the small island dynamics. 

The corresponding CZDs are bimodal, as shown in Fig. 
4b, with the first peak corresponding to small islands and 
the second peak corresponding to the large islands. This 
bimodal CZDs are completely different from those of IA 
models, where distortions of the monomodal shapes of 
ISDs and CZDs appear only for large coverages 38 . Here, 
as 9 increases, the first peak of the CZD moves to the 
left (smaller CZs for small islands) and the second one 
moves to the right (larger CZs for large islands). 

A surprising result in the inset of Fig. 4b is that the 
Gaussian right tail of the CZDs is preserved for low cov- 
erages in regime 777. Thus, in this fully reversible regime, 
the dynamics of large islands is still that of PE theory, 
which represents their competition for the capture of dif- 
fusing adatoms. The small islands are rapidly dissociat- 
ing into isolated monomers to maintain the competitive 
dynamics of the neighboring large islands. The right tails 
of the corresponding ISDs are also Gaussian, following 
the trend of compact islands already explained in Ref.— . 

In experiments, very small islands (s = 2, 3, • • •) may 
not be detected in microscopy images, except if high res- 
olution techniques are used. If that is the case, the first 
peak of the bimodal CZDs may be lost, but the Gaus- 
sian right tail will be preserved. This may suggest IA 
behavior in a system with fully reversible dynamics. 



IV. TRIANGULAR LATTICE 

Fig. 5a shows the independence of the island density 
for various values of En ■ It has the main features of the 
RA model in the square lattice, including the high tem- 
perature crossover to the fully reversible regime, where 
Ni s i and N% increase with R. However, some important 
features of RA are particular of the triangular lattice, as 
discussed below. 

The crossover from a fractal island regime I (i = 1) 
to the first regime of compact islands is equivalent to 
that in the square lattice, occurring when t% matches 
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FIG. 5: (Color online) Scaled island densities for 9 — 0.10 
and several binding energies in the triangular lattice. Dashed 
lines indicate local slopes. The insets show scaled monomer 
densities. 



T ci g and leading to the crossover relations in Eq. 
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This is illustrated in Fig. 5b. The subsequent regime is 
hereafter called 77a and corresponds to I A with i = 2, 
i. e. with the smaller stable island being a triangle with 
s = 3 adatoms, shown in Fig. 6a. From Eq. (jlip and 
the RE equations (JTJ) and ([2]) for i = 2, regime //a has 
initial scaling as 

(19) 

This trend is confirmed by the slope of the curves in Fig. 
5b after the crossover. Eq. (fl9t is consistent with a 
binding energy of the critical island E^ — E^, which is 
reasonable for a two-adatom island. 

As the temperature increases, regime Ila evolves so 
that the critical island binding energy becomes negligible 
and island and monomer densities scale as 

Nisi^RT 1 / 2 , N^R-V 2 (IIa fmal ). (20) 

This evolution parallels that observed in the square lat- 
tice. 

A crossover to another regime also occurs when the 
typical aggregation time to an island is of the order of 
the detachment time T2 of doubly-bonded atoms (Fig. 
6a). Matching Eqs. © and (JH} using fl20|), we obtain a 
crossover variable 

Y 3 = ei? 1 / 4 (21) 

and crossover scaling as 

N M ~ e 2 F 3 (y 3 ) , Ni ~ e 2 G 3 (Y 3 ) . (22) 

Fig. 5c illustrates this crossover scaling. 

In the square lattice, a similar crossover lead to regime 
III. However, here it leads to another regime with I A 
scaling due to the triangular lattice geometry. The reason 
is that an island with s = 4 atoms (Fig. 6b), is stable 
for a time much longer than the island with s = 3 atoms 
(Fig. 6a), as will be explained below. Thus a regime 
with i = 3 appears when R (or T) increases, which is 
hereafter called regime lib. Its scaling properties follow 
by matching the independence of Eqs. |T]) and © for 
i = 3 and that in Eq. ((22)) : 

N zsl ~ e- 2 / 5 iT 3 / 5 , Ni ~ e 2 / 5 iT 3 / 5 (Hb initial ). 

(23) 

The slope —12/5 indicated in Fig. 5c after the crossover 
confirms this result. Similarly to the other regimes, it 
evolves to an e-independent form as the temperature in- 
creases: 

A^,~iT 3 / 5 , N^R- 3 / 5 {IIb flnal ). (24) 

Figs. 6a and 6b show the islands with s — 3 and s = 4 
atoms loosing one adatom, which provide the (critical) is- 
lands with i = 2 and i = 3, respectively. Both processes 
of loosing one adatom occur in a characteristic time Ti 
[Eq. Q3] . The island with i = 2 needs a time ti to break 
into two monomers, while the island with i = 3 needs a 




FIG. 6: (Color online) Illustration of stable islands of size (a) 
s = 3 (i = 2) and (b) s = 4 (i — 3) on a triangular lattice 
before (left) and after (right) loosing a particle. 

time T2 to break into isolated monomers. The addition 
of these times do not change the order of magnitude of 
the mean-field rate 1/t2- Thus, in a mean- field approach 
(RE theory), only the regime with i = 2 is predicted 
in the triangular lattice. On the other hand, the sta- 
bility time of the remaining islands in Figs. 6a and 6b 
is affected if the recurrence property of two-dimensional 
random walks is considered, which drastically changes 
mean-field predictions. 

The probability that a random walker returns to the 
origin after n steps is of order 1 /n 2 — . Thus, the probabil- 
ity that the island with i = 2 breaks before the return of 
the detached atom (Fig. 6a) is of order J™ l/n 2 dn ~ e. 
Analogously, the probability that the island with i = 3 
breaks before the return of the detached atom (Fig. 6b) 
is of order J D l/n 2 dn ~ e 2 . This means that the life- 
time of the island with i — 3 is (T 1 times larger than the 
lifetime of the island with i = 2. For e <C 1, both regimes 
will be observed. 

The above probabilities suggest that the rates of ag- 
gregation and detachment are not correctly expressed in 
RE theory. Indeed, this is confirmed by simulation of 
IA models ^ 28 ' 41 . However, the success of RE theory fol- 
lows from the fact that all rates are overestimated (and 
characteristic times underestimated) by similar factors, 
so that the final scaling on R, 8 and E^ is correct. Fol- 
lowing these ideas, we incorporate the above factor e _1 
in the mean-field approach by stating that the effective 
lifetime for the critical island with i = 3 is 

r 3 - e-V 2 - 1/ (De 3 ) . (25) 

Since the monomer density decreases with temperature 
slower than island density in regime lib, it is expected to 
be followed by the fully reversible regime 777, similarly 
to the square lattice. Matching t 3 with the aggregation 
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time [Eq. ((5J] for the i = 3 regime, the crossover scaling 
variable is defined as 

F 4 = eR 2/15 . (26) 
The crossover scaling is obtained considering Eq. (j2"4")l : 

N iM ~ e 9 / 2 F 4 (y 4 ) , iVi ~ e 3 G 4 (Y 4 ) • (27) 

Fig. 5d illustrates this crossover scaling (again with some 
deviations in the collapse of N± due to corrections to 
scaling). 

After this crossover, the high-temperature, fully re- 
versible regime III is observed, with the same basic fea- 
tures of the square lattice. The island and monomer den- 
sities scale only with e as 

N lsl ~e>' 2 , iVx-e 3 {III), (28) 

thus increasing with temperature. This is also a steady 
state, nearly deposition independent, where N± follows 
from matching of aggregation and detachment times 

[Eqs. (|8]) and ([25])] and N isl ~ N^ 2 , a relation char- 
acteristic of the largest available critical island. 

The features of ISDs and CZDs in the triangular lat- 
tice are similar to those of the square lattices, including 
the Gaussian tails observed at high temperatures and ex- 
tending to the fully reversible regime. 

V. CONCLUSION 

We analyzed the scaling of island density, CZDs, and 
ISDs in a model of fully reversible island growth where 
the activation energy for an adatom hop depends only 
on the number of NNs at its initial position (Clarke- 
Vvedensky model). A scaling approach based on rate- 
equation results for IA models was presented and sup- 
ported by numerical simulations in square and triangular 
lattices. 

At low temperatures, a regime / with fractal islands is 
observed, corresponding to I A with critical island size 
i = 1. CZDs and ISDs are usually intermediate be- 
tween those of RSA (negligible diffusion of free adatoms) 
and those of PE theory (competition of large islands for 
rapidly moving adatoms), except for large values of the 
binding energy En, where the regime is extended up to 
R ~ 10 7 and a nearly Gaussian right tail of CZDs ap- 
pears. 

A crossover to a second regime is always observed 
as the reduced model parameters satisfy ei? 2 / 3 ~ 1, in 



agreement with previous works. In the square (trian- 
gular) lattice, a regime with scaling similar to IA with 
i = 3 (i = 2) is observed after the crossover. The ini- 
tial dependence of island density and monomer density 
with e disappears as the temperature increases, so that 
those quantities scale only with R close to the subsequent 
crossover. 

In the triangular lattice, there is a crossover to a third 
I A regime, with i — 3. This is explained by the recur- 
rence properties of random walks in two dimensions, since 
the detachment times of doubly bonded atoms in islands 
with 3 and 4 atoms are the same and that feature could 
not be predicted solely with a mean-field approach. 

Although the CZDs shows some deviations from the 
WS in these regimes with i > 1, their right tails are 
typically Gaussian, indicating that the PE theory is rea- 
sonable to represent the competition of large islands for 
free adatom capture. Exceptions appear for small values 
of En, where those regimes may appear for small R. 

A final crossover to a fully reversible regime III is ob- 
served as Y 2 ~ 1 (Y 4 ~ 1) in the square (triangular) 
lattice. In contrast to IA scaling, this regime is charac- 
terized by an increase of the island density and of the 
monomer density with temperature, for fixed coverage. 
A large density of small islands coexists with a small 
density of very large islands. The corresponding ISDs 
show a high left tail and the CZDs are bimodal. The 
Gaussian right tails of CZDs and ISDs are also observed, 
despite the reversible nature of the aggregation, and are 
explained because the small islands do not change the 
main mechanisms for the large islands to compete for 
free adatom aggregation. 

These results show that some predictions of PE theory 
for the CZDs are kept in reversible aggregation, and so 
are some consequences for the ISDs (e. g. the Gaussian 
right tails) . In experiments where the resolution of imag- 
ing techniques is not high, detection of small islands may 
be difficult and even the high temperature regimes may 
have ISDs and CZDs similar to those of IA. This may ex- 
plain the success of PE predictions to describe a variety 
of experimental results, for various materials and growth 
conditions^— . 
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